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ABSTRACT 

The next generation of galaxy surveys will attempt to measure the baryon oscilla- 
tions in the clustering power spectrum with high accuracy. These oscillations encode 
a preferred scale which may be used as a standard ruler to constrain cosmological 
parameters and dark energy models. In this paper we present simple analytical fit- 
ting formulae for the accuracy with which the preferred scale may be determined in 
the tangential and radial directions by future spectroscopic and photometric galaxy 
redshift surveys. We express these accuracies as a function of survey parameters such 
as the central redshift, volume, galaxy number density and (where applicable) pho- 
tometric redshift error. These fitting formulae should greatly increase the efficiency 
of optimizing future surveys, which requires analysis of a potentially vast number of 
survey configurations and cosmological models. The formulae are calibrated using a 
grid of Monte Carlo simulations, which are analyzed by dividing out the overall shape 
of the power spectrum before fitting a simple decaying sinusoid to the oscillations. 
The fitting formulae reproduce the simulation results with a fractional scatter of 7% 
(10%) in the tangential (radial) directions over a wide range of input parameters. We 
also indicate how sparse-sampling strategies may enhance the effective survey area if 
the sampling scale is much smaller than the projected baryon oscillation scale. 
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1 INTRODUCTION 

Baryon oscillations in the galaxy power spectrum have re- 
cently emerged as a promising standard ruler for cosmol- 
ogy, potentially enabling precise measurements of the dark 
energy parameters with a minimum of systematic errors 
(Cooray et al. 2001; Eisenstein 2002; Blake & Glazebrook 
2003; Seo & Eisenstein 2003; Linder 2003). The large- 
scale linear clustering pattern contains a series of small- 
amplitude, roughly sinusoidal, modulations in power of iden- 
tical physical origin to the acoustic peaks observed in the 
Cosmic Microwave Background (CMB) (see Eisenstein & 
Hu 1998; Meiksin, White & Peacock 1999 and references 
therein). These oscillations encode a characteristic scale - 
the sound horizon at recombination - which can be accu- 
rately calibrated using the linear physics of the CMB. The 
apparent value of this preferred scale, deduced from a slice 
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of a galaxy spectroscopic redshift survey, depends on the 
assumed cosmological distance and expansion rate at the 
slice redshift z, which control the mapping of redshifts to 
physical co-ordinates in the tangential and radial directions, 
respectively. The baryon oscillations can therefore be used 
to measure the angular diameter distance Da{z) and Hub- 
ble parameter H{z) in units of the sound horizon, over a 
series of redshift slices. The acoustic signature may also be 
measured from photometric redshift surveys (Blake & Bridle 
2005), although the smearing of radial information implies 
that only Da{z) may be determined with any confidence. 
The preferred scale was recently identified in the clustering 
pattern of Luminous Red Galaxies in the Sloan Digital Sky 
Survey (Eisenstein et al. 2005) and used to constrain cos- 
mological parameters for the first time. A power spectrum 
analysis of the final 2-degree Field Galaxy Redshift Survey 
produced consistent results (Cole et al. 2005). 

A number of techniques have been employed to estimate 
the accuracy with which the baryon oscillation scales may be 
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determined by future galaxy surveys. Blake & Glazebrook 
(2003, hereafter BG03) used a Monte Carlo, semi-empiricaJ 
approach in which realizations of spectroscopic redshift sur- 
veys were created from an underlying linear power spec- 
trum. The acoustic scale weis recovered for each realization 
by first dividing out the overall shape of the measured power 
spectrum, then fitting a simple empirically-motivated decay- 
ing sinusoid to the baryon oscillations, up to a maximum 
wavenumber determined by a conservative estimate of the 
extent of the linear regime at the redshift in question. The 
scatter in the best-fitting values of the acoustic 'wavelength' 
across the realizations represents the accuracy with which 
the preferred scale may be extracted in such an experiment. 
Blake & Bridle (2005) extended this methodology to photo- 
metric redshift surveys. 

A feature of this 'model-independent' approach is that 
the information contained in the oscillations is decoupled 
from that encoded by the overall shape of the power spec- 
trum, which is divided out prior to fitting the sinusoid, and 
which may be subject to smooth broad-band systematic tilts 
from such effects as poorly-modelled redshift-space distor- 
tions, scale-dependent bias and non-linear growth of struc- 
ture. On the other hand, the power spectrum shape does 
also depend on the cosmological parameters, and combined 
measurements in the tangential and radial directions permit 
an Alcock-Paczynski test (Yamamoto, Bassett & Nishioka 
2005). BG03 discard this potentially useful information for 
the benefit of simulated measurements that axe more 'ro- 
bust' against the presence of systematic errors. Furthermore 
the BG03 analysis contains various approximations, as dis- 
cussed below. 

Various other studies predicting the standard ruler ac- 
curacies from baryon oscillation surveys have fitted a full 
power spectrum template in order to estimate the cosmic 
distance scale. Several papers (Sco & Eisoustcin 2003; Hu 
& Haiman 2004; Amendola, QuerccUini & Giallongo 2005; 
Huetsi 2005) employ Fisher matrix techniques to recover 
predicted errors in the cosmological quantities. In addition, 
galaxy catalogues extracted from N-body simulations, and 
therefore incorporating realistic non-linear and galaxy bi- 
asing efi^ects, have been analyzed by Angulo et al. (2005), 
Springel et al. (2005), Seo & Eisenstein (2005) and White 
(2005). In these studies the power spectrum is typically 
'linearized' using a polynomial function, and is then fitted 
with a linear-regime model using normal clii-squared tech- 
niques. Whilst the initial results from these investigations 
are encouraging, it appears that redshift-space distortions 
slightly degrade the baryon oscillation accuracy in the ra- 
dial direction when fitting a full power spectrum template 
(Seo & Eisenstein 2005). Larger simulations are required 
in order to quantify accurately the potential infiuence of 
low-level systematic errors from redshift-space distortions, 
scale-dependent bias and non-linear growth of structure in 
the baryon oscillations technique. 

The 'model-independent' and 'full template' methods 
are complementary, with the minimal sinusoid-fitting pro- 
viding an efi'ective lower limit to the efficacy of the tech- 
nique, and the full power spectrum shape fit indicating what 
may be achieved with more assumptions. In fact, in the 
regime where the oscillations are measured with high star 
tistical confidence, they encode most of the potential for 
constraining cosmology (Hu & Haiman 2004) and the accu- 



racies predicted by the two techniques agree reasonably well 
(Glazebrook & Blake 2005). 

The simplicity of the 'model-independent' technique of 
BG03 implies that the resulting baryon oscillation accuracies 
scale in a relatively predictable manner with the spectro- 
scopic and photometric survey parameters: central redshift, 
survey volume, galaxy number density and (where applica- 
ble) photometric redshift error. The purpose of this study is 
to provide accurate fitting formulae for these standard ruler 
accuracies in terms of the survey parameters. We consider 
baryon oscillation measurements in both the tangential and 
radial directions for spectroscopic redshift surveys, and in 
the tangential direction alone for photometric redshift sur- 
veys. These formulae will be considerably more efficient to 
implement than full Monte Caxlo power spectrum realizar 
tions, and should prove useful for planning future galaxy 
surveys with the goal of measuring the dark energy param- 
eters. The general design optimization of such surveys in- 
volves consideration of a potentially vast parameter space 
of survey configurations and cosmological models (see Bas- 
sett 2005; Bassett, Parkinson & Nichol 2005), and in this 
context a fitting formula is invaluable. 

We assume a fiducial ACDM fiat cosmological model 
with matter density fim = 0.3 although, as discussed in 
Section 5, our results apply more generally. 



2 SIMULATED GALAXY SURVEYS 
2.1 Starting assumptions 

We begin by emphasizing some approximations inherent in 

our analysis: 

• We utilize no information encoded in the overall shape 
of the power spectrum, which is divided out by a smooth 
polynomial prior to the baryon oscillation fit. 

• We assume that the power spectrum errors can be de- 
scribed by (correlated) Gaussian statistics for wavenumbers 
fc up to a maximum fcmax, specified by a conservative esti- 
mate of the extent of the linear regime at the redshift in 
question, and that modes with scales k > fcmax provide no 
information. 

• We employ an approximate parameterized fit (a decay- 
ing sinusoid) for the baryon oscillation signature (sec BG03) . 
We neglect the small scale-dependent phase shifts of the 
acoustic peaks and troughs. 

• We assume that shot noise can be described by Poisson 
statistics at the per cent level. This may not be the case, and 
planned galaxy surveys may require slightly higher galaxy 
number densities in order to be cosmic-variance limited. 

The measurements of future galaxy surveys should be 
fitted with accurate templates, marginalizing over model un- 
certainities, rather than empirical sinusoids. However, we 
believe that the approximations contained in our current 
analysis are acceptable because: 

• Comparing our results with full Fisher matrix simula- 
tions, such as those of Seo & Eisenstein (2003, 2005), our 
inferred standard ruler accuracies are comparable (albeit 
30 — 50% larger, refiecting our more conservative approach). 

• A major simulation effort is still required to model halo 
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bias, non-linear structure formation and redshift-space dis- 
tortions to the required accuracy (together with their effects 
on power spectrum mode correlations). Using this informal 
tion in our templates at this stage would create the possi- 
bility of additional systematic error. 

• Our approach enables us to explore a very large number 
of survey configurations and cosmological models. 



2.2 The simulation grid 

In order to explore the scalings of baryon oscillation accura- 
cies with galaxy survey parameters, we created a largo grid 
of simulated spectroscopic and photometric rodshift surveys. 
For each survey configuration we used Monte Carlo realiza- 
tions to determine the accuracy with which the standard 
ruler could be measured, as described below. For spectro- 
scopic surveys, wc derived accuracies in the tangential and 
radial directions. For photometric surveys, the damping of 
radial information implies that baryon oscillations may only 
be measured in the tangential direction (Seo & Eisenstein 
2003, Blake & Bridle 2005). 

In cosmological terms, the measurement accuracies of 
the tangential and radial baryon oscillation scales determine 
the precision with which the quantities r{z)/s and r'{z)/s 
may be inferred, where r(z) is the co-moving distance to 
the redshift slice z, r'(z) = dr/dz — c/H{z) where c is the 
speed of light and H{z) is the Hubble parameter measured 
by an observer at redshift z, and s is the (co-moving) sound 
horizon at recombination. Given that s = 6a "X r(CMB), 
where 9a is the (accurately-known) angular scale of the first 
CMB acoustic peak and r(CMB) is the inferred distance to 
the surface of last scattering, then the acoustic oscillations 
may be thought of as measuring the quantities r{z)/r{CMB) 
and r'(2)/r(CMB) (see Eisenstein et al. 2005). 

Our grid of Monte Carlo simulations was generated by 
varying four survey parameters: the central redshift z, the 
survey area A (in 10^ deg^), the survey width Sz (such that 
the survey ranges between redshifts z — 5z and z + 6z) and 
the number density of galaxies n (in 10"'' ft'' Mpc^'^). For 
photometric redshift surveys, a fifth parameter was added; 
the r.m.s. error in redshift cr^, expressed via the parameter 
fo = Cz/{i + z). Given the potential complexity of target 
selection techniques for these surveys, we do not consider re- 
alistic galaxy redshift distributions, but instead populate the 
survey volume uniformly with number density n. This will 
be a good approximation for any relatively narrow survey 
redshift slice (or for high-enough number density the mear- 
surements will be cosmic-variance limited and independent 
of n) . Further details of our simulations of spectroscopic and 
photometric redshift surveys are given below. 



2.3 Spectroscopic redshift surveys 

The trial parameter values for the simulated spectroscopic 
redshift surveys were: 

z = 0.2 -> 3.4 in steps of 0.2 

A = 0.1,0.25,0.5,0.75,1.0,1.5,2.0,3.0,4.0,5.0 

Sz = 0.025,0.05,0.075,0.1,0.125,0.15,0.2,0.3,0.4,0.5 

n = 0.04,0.08,0.16,0.32,0.64,1.3,2.6,5.1,10.2,20.6 



i.e., 16,600 configurations were analyzed (given that cases 
with 6z > z are excluded). 

The Monte Carlo methods utilized to analyze the simu- 
lated spectroscopic redshift surveys were similar to those em- 
ployed by BG03 and described in more detail by Glazebrook 
& Blake (2005). However, some improvements in speed were 
required to process such a large grid of surveys. In BG03, 
Monte Carlo realizations of surveys were generated by per- 
forming Gaussian realizations of the underlying power spec- 
trum, and Poisson sampling the resulting density fields. The 
covariance matrix of the power spectrum bins was effec- 
tively determined numerically by averaging over the Monte 
Carlo realizations. This process is too time-consuming for 
exploring a large grid of surveys. Therefore in this study 
we determined the covariance matrix analytically from the 
Fourier transform of the survey window function, by evalu- 
ating the sums given in Feldman, Kaiser & Peacock (1994), 
who present an optimal estimator for the power spectrum. 
Our analysis assumes a full conical survey geometry, and 
hence includes the convolution of the underlying power spec- 
trum with the survey geometry and the correlations between 
the Fourier bins. 

In more detail: for a given survey configuration we en- 
closed the survey cone (sampled by a uniform number den- 
sity n) with a cuboid of volume Vcub, and determined the 
Fourier transform of the window function, which we write 
as Wk- We assumed a model linear power spectrum given 
by the fitting formula of Eisenstein & Hu (1998), using fidu- 
cial cosmological parameters Qm = 0.3, Qa = 0.7, h = 0.7, 
ilb/^m = 0.15, Us — 1 and as = 1. We scaled this power 
spectrum to redshift z using the linear growth factor D{z) 
of Carroll, Press & Turner (1992). We convolved the in- 
put power spectrum with the survey window function. We 
restricted our analysis to scales larger than a maximum 
wavenumber k corresponding to a conservative estimate of 
the transition scale between the linear and non-linear clus- 
tering regimes (see BG03) . Our surveys assume a linear bias 
factor bo = 1 for galaxies with respect to matter, which is 
likely to be conservative at high redshifts. Our fitting for- 
mulae can be simply adapted for 6o 7^ 1 as explained below. 

In order to obtain the power spectrum covariance ma- 
trix, we note that equation 2.5.2 in Feldman, Kaiser & Pea^ 
cock (1994) reduces to the expression (see also Tadros & 
Efstathiou 1996): 

_ < SPjSPj > _ 2 J2k,k' IW^fc-fc'l" 

where P = P{k) is the power spectrum amplitude typical of 
bins i and j, the summation in the numerator is evaluated 
between all A'sum separate pairs of modes (fe, k') in bins i and 
j, and the summation in the denominator is evaluated over 
all Fourier modes and is equal to (V/Vcub)^ for a window 
function which is either a uniform value or zero. The factor 
of 2 reflects the fact that only half of the measured Fourier 
modes are independent, owing to the reality condition of the 
density field. For a survey in a uniform box, Wk = unless 
fc = and equation 1 reduces to Ca — 2/m, where rn is the 
total number of Fourier modes in bin i. We evaluated the 
double sum in equation 1 using a Monte Carlo integration 
scheme. We binned the power spectrum into a 2-dimensional 
grid of radial and tangential components (if the x-axis is the 
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radial direction, then ferad = and fctan = \/ky-\- fci). We 
used Fourier bin widths Afcrad = Afctan = 0.01 /i Mpc~^, 
unless a survey dimension L was sufficiently small that the 

corresponding spacing of the Fourier modes 2ti / L > 0.005 h 
Mpc~^, in which case we set A/c = A-k/ L (i.e., the minimum 
thickness of a bin in our analysis is 2 Fourier modes) . Given 
that the acoustic 'wavelength' in Fourier space is fcA ~ 0.06 h 
Mpc~^, we rejected a survey configuration if A/c > 0.03 /i 
Mpc~^. We tested our code using an analytically-tractable 
survey window function for which the sums in equation 1 
could be evaluated in closed form. 

Having determined the covariance matrix, we created 
many Gaussian realizations of correlated power spectrum 
measurements using the technique of Cholesky decomposi- 
tion. The acoustic 'wavelengths' in the tangential and radial 
directions were fit to these realizations via a simple empiri- 
cal decaying sinusoid, using the same method as Glazebrook 
& Blake (2005), and the accuracies with which these scales 
could be measured was inferred using the scatter in the 
resulting best-fitting 'wavelengths' across the realizations. 
Tests repeating the analysis for identical survey configura- 
tions showed that the scatter in determination of the stan- 
dard ruler accuracy y owing to numerical noise, due to the 
approximate summation scheme for equation 1 and the finite 
number (400) of power spectrum realizations, was about 5% 
of y (7% of y) in the tangential (radial) direction. 



fitting formula of Eisenstein & Hu (1998) , and m is the num- 
ber of Fourier modes contributing (i.e. those contained in 
an annulus of radius fctan, radial thickness 5fctan and length 
<5fcrad). The value of m was determined from the density-of- 
states in fe-space, pk = F/(27r)^: m = Pk x 27rfetan x 5fctan x 

For each survey configuration, we generated a large 
number of Monte Carlo realizations of noisy power spec- 
tra by adding a Gaussian variable of standard deviation 5P 
(given by equation 2) to the damped model power spectrum 
Pexp [— (feradCTr)^]. We binned each power spectrum real- 
ization in tangential Fourier bins of width (5fctan ~ 0.01 h 
Mpc^^ by averaging cells in the radial direction up to a 
maximum value of fcrad = l-B/cTr, beyond which the power 
spectrum contains very little signal owing to the damping. 
For each realization, we then divided the binned power spec- 
trum by a smooth 'reference spectrum' and fitted the result 
with the empirical decaying sinusoid used in BG03. As in the 
case of the spectroscopic redshift surveys, the scatter in the 
best-fitting values of the sinusoidal 'wavelength' across the 
Monte Carlo realizations was taken as the tangential baryon 
oscillation accuracy for this survey configuration. Tests re- 
peating the analysis for identical configurations showed that 
the scatter in the determination of the standard ruler accu- 
racy y owing to numerical noise was about 5% of y. 



2.4 Photometric redshift surveys 



3 THE FITTING FORMULA 



The trial parameter values for the simulated photometric 
redshift surveys were: 

z = 0.5,0.75,1.0,1.25,1.5,1.75,2.0,2.5,3.0 

A = 0.5,0.75,1.0,2.0,5.0,10.0 

Sz = 0.1,0.2,0.3,0.4,0.5 

n = 0.1,0.2,0.5,1.0,2.0,5.0,10.0 



0.01,0.02,0.03,0.04, 0.05 



i.e., 9,450 configurations were analyzed. 

For this set of Monte Carlo simulations a flat-sky ap- 
proximation was assumed, as described by Blake & Bridle 
(2005). Therefore, unlike for the spectroscopic survey anal- 
ysis, no window function effects are considered. The photo- 
metric redshift error distribution was assumed to be a Gaus- 
sian function, hence the power spectrum was assumed to be 
damped in the radial direction by a factor eXp [ (fcradC^r) ], 
where (Jr = ao{l+z) xdr/dz (Blake & Bridle 2005, equations 
3 and 4). 

In order to speed up the computation in the same style 
as described above for spectroscopic surveys, the error in the 
power spectrum measurement for photometric redshift sur- 
veys was not determined via radially-smeared Monte Carlo 

realizations of the density field (as in Blake & Bridle 2005) 
but instead using an analytical approximation. The error 
in the power spectrum measured in a Fourier cell centred 

at (fctan, fcrad) of width (5fctan , 5fcrad ) was assumcd to be the 
usual combination of cosmic variance and shot noise: 



1 / 2 1 

SP = —= Pexp [- (fcrad O-r) ] + " 



(2) 



where P = P{k) = P{y/k^^^ + fc^ad) the (undamped) 
value of the model power spectrum in the bin, taken from the 



In this Section we develop an analytic expression for the 
accuracy of measurement of the baryon oscillation scale in 
terms of the survey configuration: central redshift z, total 
volume V (in Gpc^), average number density of galax- 
ies n (in 10~^ Mpc~^) and, in the case of photometric 
redshift surveys, the r.m.s. error in co-moving co-ordinate 
tJr (in Mpc). The fitting formula contains free parame- 
ters whose values are calibrated using the grids of simulated 
surveys described in Section 2. Different fitting formula co- 
efficients were derived for standard ruler accuracies result- 
ing from spectroscopic redshift surveys (separately for the 
tangential and radial directions) and photometric redshift 
surveys (in the tangential direction only) . 

As a first approximation for the fitting formula, we as- 
sumed that the accuracy x with which the acoustic scale can 
be measured is proportional to the average fractional error 
5P/P in the power spectrum, given by the usual sampling 
formula (e.g. Tegmark 1997): 



6P 



1 



1 



(3) 



where m is the total number of independent Fourier modes 
contributing to the measurement, and P = P{k*) (in 
Mpc"*) is the value of the power spectrum amplitude at an 
average scale fc* ~ 0.2 h Mpc^^ characteristic of the baryon 
oscillations. The two terms in equation 3 represent the ef- 
fects of cosmic variance and shot noise, respectively. Given 
a fixed amount of observing time, optimal measurements of 
the power spectrum follow from a survey of depth such that 
nP ~ 1. This requirement is readily achieved by ~ 1 hr 
integrations with 8-m class ground-based telescopes. 

The number of measured Fourier modes scales with the 
total survey volume V, which determines the density-of- 
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Table 1. Best-fitting coefiicients for the fitting formula, defined by equations 6 to 9, for the three types of standard ruler accura<;y. The 
first set of parameters, from xq to b, apply only to high-accuracy baryon oscillation measurements where the simple scalings described 
by equation 6 are valid (i.e., x <S xt). The second set of parameters, from p to /3, describe the worsening standard ruler performance in 
the regime where the oscillations are just being resolved (see equation 8). The r.m.s. difference in the predictions of the fitting formulae 
and Monte Carlo simulations is listed separately for measurements of all accuracies {y < 10%) and of just high accuracy {y < 2%). 



Parameter 




Spec-2 


Spec-z 


Photo-z 






Tangential 


Radial 


Tangential 


xo (per cent) 




0.85 


1.48 


1.23 


no(xlO-3 Mpc-3) 




0.82 


0.82 


0.71 






1.4 


1.4 


1.4 


7 




0.5 


0.5 


0.61 


6 




0.52 


0.52 


0.52 


P 




2 


2 


4 


a 




7.3 


10.6 


4.2 


a 




0.26 


0.49 


0.11 


P 




0.27 


1.00 


0.42 


r.m.s. error in fitting formula {% 


< 10%) 


7.1% of J/ 


9.9% of y 


6.9% of y 


r.m.s. error in fitting formula ( 


y < 2%) 


5.5% of y 


7.4% of y 


4.0% of y 



states in Fourier space pk- m oc. pk (x V. In addition, m 
is proportional to the contributing volume in fe-space. For a 
fiducial survey spanning 1000 dcg^ from 2 = 0.6 to ^ = 1.4 
(i.e. A = 1, 6z ^ 0.4), we find V = 2.16 /i"^ Gpc=^ for 
our fiducial cosmological parameters. The number of inde- 
pendent modes contained in an angle-averaged Fourier bin 
at fe = 0.2 ft Mpc"^ of width 6k = 0.01 /i Mpc"^ is then 
m « 2.2 X lO'*, yielding a measurement of the power spec- 
trum in this bin with an accuracy of about 1% using equa- 
tion 3. As the baryon oscillations have a fractional amplitude 
of roughly 5%, this constitutes a high- significance detection. 

For photometric redshift surveys, the radial smearing 
damps out the useful signal in all Fourier modes with small- 
scale radial fe- values ferad ^ l/o"r. Therefore, m oc l/ur. For 
a typical photometric redshift performance, l/fr,. ~ 0.02 ft 
Mpc~^ . Comparing this scale with the extent of the available 
linear regime, fe<0.2ft Mpc^^, we find that the number of 
usable Fourier modes is diminished by roughly an order of 
magnitude for photometric surveys. 

If we also include the scaling of the power spectrum 
with redshift as the linear growth factor D{z) and a linear 
bias parameter feo, we can re- write equation 3 as: 



1 + 



Wcff -P(zo)^ 



(4) 



where n^ff (~ 1/P(fc*) at z = zo) is a fiducial number density 
(a fitted parameter) and D{zo) is the linear growth factor 
at a fiducial redshift z = zo- The factor y/a^ only appears 
for photometric redshift surveys. The assumption of a linear 
scalc-indcpcndont bias factor 60 will be incorrect in detail. 
However, given that in our simulations the overall shape of 
the power spectrum is divided out prior to fitting the sinu- 
soidal function, our results are not sensitive to such details. 
The value of 60 may be interpreted as the boost in the power 
spectrum of galaxies with respect to that of dark matter at 
the characteristi c scale of the bary on oscillations, fe* ~ 0.2 ft 
Mpc~^: bo = ^/ Pga.i{k*) / Pdmik*). We note that the quan- 
tity 60 X D{z) in equation 4 is observed to be roughly con- 
stant with redshift for L* galaxies (Lahav et al. 2002). 



We normalize equation 4 to an accuracy x = xo for a 
fiducial survey of volume V = Vo and (where applicable) 
redshift error ar = cr,. n: 



X = Xo 




(5) 



We take zo = 1 and the same fiducial survey as above 
{A = 1, 5z = 0.4) such that Vo = 2.16 ft"^ Gpc^ and 
D{zo) = 0.61 (Carroll, Press & Turner 1992). We also as- 
sume o'j.,() — 34.1 h^^ Mpc, corresponding to a redshift error 
az — 0.01(1 + z) at z = zo- 

The fitting formula must also reflect the increase with 
redshift of the extent of the linear regime (i.e. the number of 
acoustic peaks which may be fitted), which enables a more 
accurate determination of the acoustic scale for a fixed power 
spectrum precision (see Glazebrook & Blake 2005, Figure 7). 
This is accomplished via an empirical power-law in z, which 
is cut off at a maximum redshift z = Zm a-t which all of the 
high-amplitude peaks are visible: 



= Xo 




n hlD{zY 
11 blD{zY 



Z < Z-a 



Z> Zjn (6) 



where 7 > is a fitted parameter. In addition, given that 
the amplitude of the power spectrum decreases with increas- 
ing fe, the variation in the extent of the linear regime with z 
changes the average amplitude of P{k) included in the anal- 
ysis, and hence the value of rieff . We described this variation 
by 



rioff = no 
= m 



1-6 



\ ZjnJ 



Z < Zu 
Z> Zn 



(7) 



The fitting formulae of equations 6 and 7, containing five free 
parameters {xo, no, b, 7, Zjn), work well for high-precision 
measurements of the acoustic scale. However, in the regime 
where the oscillations are just being resolved, the scaling of 
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the 'model-independent' accuracy with (for example) survey 
volume is more rapid than V~^^'^. We therefore modified the 
accuracy a; to a new value y where 



where Xt is a characteristic accuracy and p > is a free 
parameter. Equation 8 is designed such that y ^ x as, x ^ Q. 
The quantity y is the final predicted standard ruler accuracy 
of the fitting formula. 

Empirically, we found that the quantity xt has a depen- 
dence on survey volume V and redshift z: 

^* = '^(^) (t)' 

f V \" 

= a (^^j 2 > 2m (9) 

Equations 8 and 9 hence introduce another four parameters 
(p, a, a and /3) which describe the departure from the simple 
scaling in the high-accuracy regime. 

We emphasize that planned baryon oscillation surveys 
should aim to reach the regime in which the acoustic features 
have been properly resolved and the high-accuracy scalings 
of equation 6 apply (i.e., x <C xt). Otherwise, the detection 
of oscillations will be of poor significance in many Monte 
Carlo realizations of the planned survey. The modifications 
represented by equations 8 and 9 are included to ensure that 
equation 6 is not applied in the regime where the oscillations 
are poorly detected, which would result in over-optimistic 
predictions of the standard ruler accuracies. Equation 6 has 
some elements in common with the formula suggested by 
Bernstein (2005, equation 42) for the tangential baryon os- 
cillation accuracy. 



4 FITTING FORMULA COEFFICIENTS 

We varied the free parameters of the fitting formulae to ob- 
tain the best fit to the grid of Monte Carlo simulated sur- 
veys (in the sense of the lowest standard deviation of the 
fractional variation). If the Monte Carlo accuracy of the 
acoustic scale was poorer than 10% for a grid point, then 
that survey configuration was assumed to provide no mea- 
surement of the baryon oscillations and was ignored in the 
fitting process. Of the 16,600 spectroscopic survey config- 
urations, 5,441 (4,576) were included in the determination 
of the fitting formula coefficients in the tangential (radial) 
direction. Of the 9,450 photometric survey configurations, 
3,765 were included. 

Results are listed in Table 1 for the cases of spec- 
troscopic surveys (tangential and radial directions) and 
photometric surveys (tangential direction). Considering the 
whole regime of standard ruler accuracies better than 10% 
(y < 0.1), the r.m.s. difference in the predictions of the fit- 
ting formulae and Monte Carlo simulations is about 7% 
of y (10% of y) in the tangential (radial) direction. For 
high-precision measurements with accuracies better than 2% 
{y ~ X < 0.02) the formulae perform significantly better: 
in this regime the fitting formulae recover the baryon os- 
cillation accuracies to better than ±0.1%. Given that the 
numerical noise in the grid of simulated surveys, resulting 
from the Monte Carlo realizations, is approximately 5% of 



y (7% of y) in the tangential (radial) directions, thus con- 
stituting a significant fraction of the scatter, these fitting 
formulae perform remarkably well. 

The parameter Zm, which is the redshift at which the 
improvement in the baryon oscillation accuracy (for fixed 
survey volume and number density) saturates, was con- 
strained to have the same value for all types of survey. For 
the best-fitting value, Zm = 1.4, a conservative estimate of 
the extent of the linear regime (see BG03) is fcun = 0.25 h 
Mpc~^, encompassing essentially the whole range of high- 
amplitude acoustic peaks. The parameters no and b were 
also constrained to be equal in the tangential and radial 
directions for spectroscopic surveys. The best-fitting value, 
no = 8.2 X 10"'' Mpc~^, corresponds to an effective power 
spectrum P --^ 1/no « 1200 ft^^ Mpc"'. This is very reason- 
able, given that for z = zo — 1 (see equation 4) the am- 
plitude of the power spectrum is P{k, z) = P{k, 0)D{z)^ Ri 
1200 h-^ Mpc^ at fe = 0.19 h Mpc'^. 

Comparing the accuracies of measuring the tangential 
acoustic scale with spectroscopic and photometric galaxy 
surveys we find that, assuming an identical number density 
and redshift range, the photometric survey (with r.m.s. er- 
ror in radial co-ordinate ar) must cover an area exceeding 
the spectroscopic survey by a factor ~ 2.1(crr/34.1 Mpc) 
to produce the same level of tangential accuracy. Blake & 
Bridle (2005) present a wider range of comparisons. 

Figures 1 to 3 compare the accuracies predicted by 
the fitting formulae with those obtained from the grid of 
simulated surveys for the three types of baryon oscillation 
measurement, illustrating the tightness of the fits. Figure 
4 plots a histogram of the fractional difference in the fit- 
ting formula and Monte Carlo accuracies for the case of the 
tangential acoustic scale from spectroscopic surveys, demon- 
strating that the scatter approximately follows a Gaussian 
distribution. There is a small systematic offset in the mean 
difference; when determining the fitting formula coefficients 
we require that this offset is less than 3% of y. 

In order to demonstrate further the performance of the 
fitting formulae, Figure 5 compares the predictions of the 
formulae with the Monte Carlo data points for measure- 
ments of tangential baryon oscillations from spectroscopic 
surveys at « = 1, as a function of survey volume. The vari- 
ous curves (and point styles) correspond to different values 
of number density n. The agreement in the shape and offset 
of the curves is excellent. 

Figure 6 plots the fitting formulae accuracies for spec- 
troscopic surveys against pairs of survey parameters: {A, n) 
and {z,Sz). Figures 7 to 10 plot the accuracies against sur- 
vey area for some more specific configurations of interest. For 
these last four figures we assume a number density of objects 
such that shot noise is unimportant ('nP = 3', see Glaze- 
brook & Blake 2005 Figure 1 for the required number density 
as a function of redshift). Figure 7 considers spectroscopic 
redshift surveys covering the redshift ranges 0.5 < « < 1.3 
and 2.5 < z < 3.5, which are naturally probed using op- 
tical spectrographs. Figures 8 and 9 display tangential and 
radial accuracies for a more general range of spectroscopic 
survey configurations in redshift slices of thickness 0.2 from 
z = 0.4 to z = 2.0. As redshift increases, the gain in survey 
volume with z saturates and thus the curves converge. Fig- 
ure 10 plots baryon oscillation accuracies from photometric 
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Figure 6. Dependence of the fitting formulae accuracies for spectroscopic surveys on pairs of survey parameters. For the left-hand plots 
we vary z and Sz, fixing A = 2 and n = 5.1. For the right-hand plots we vary A and n, fixing z = I and Sz = 0.5. The accuracy is plotted 
as l/y. 
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Figure 1. Comparison of the fitting formula and Monte Carlo 
simulation accuracies of measuring the tangential acoustic scale 
from spectroscopic redshift surveys. A significant fraction of the 
scatter is due to numerical noise in the simulations. 




Fitting formula accuracy (%) 

Figure 2. Comparison of the fitting formula and Monte Carlo 
simulation accuracies of measuring the radial acoustic scale from 
spectroscopic redshift surveys. A significant fraction of the scatter 
is due to numerical noise in the simulations. 



redshift surveys (with redshift error parameter ao = 0.03) 
in redshift slices of thickness 0.5 from z — 0.5 to 2; = 3.5. 

Comparing the predictions of the fitting formulae with 
results from the full Monte Carlo method of BG03 (e.g. 
Glazebrook & Blake 2005, Table 1) we find that the mean 
difference is about 5% of y and the standard deviation of the 
difference is roughly 10% of y. We can also compare the fit- 
ting formulae prediction with the accuracy of measurement 



of the acoustic scale by the SDSS LRG sample (Eisenstein 
et al. 2005). This survey covers sky area 3816 deg^ and red- 
shift range 0.16 < z < 0.47 {V = 0.72 h'^ Gpc^). The 
galaxy number density varies with redshift, but we take an 
effective value UcS = 10~* Mpc~^ and a galaxy bias cor- 
responding to (78 — 1.8 (Eisenstein et al. 2005). The fitting 
formulae predict measurement accuracies of 6.4% (8.5%) in 
the tangential (radial) direction, using just the oscillatory in- 
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Figure 3. Comparison of the fitting formula and Monte Carlo 
simulation accuracies of measuring the tangential acoustic scale 
from photometric redshift surveys. A significant fraction of the 
scatter is due to numerical noise in the simulations. 
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Figure 5. Comparison of the fitting formula (lines) and Monte 
Carlo simulation accuracies (data points) of measuring the tan- 
gential acoustic scale from spectroscopic surveys at 2 = 1, plotted 
as a fuuc:liou of survey volume. The different lines and data point 
styles correspond to the 10 different values of the number density 
n listed in Section 2.3, ranging from low density (upper right) 
to high density (lower left). The agreement between the fitting 
formula and Monte Carlo simulations is excellent. 
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Figure 4. Distribution of the fractional difference (Ay/y) be- 
tween the fitting formula and Monte Carlo simulation accuracies 
of measuring the tangential acoustic scale from spectroscopic red- 
shift surveys. The scatter is well-described by a Gaussian distri- 
bution, with very few outliers. The overall r.m.s. difference in 
accuracies is 7% of y (see Figure 1) . There is a small offset in the 
mean difference (2% of y). 
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Figure 7. Tangential and radial baryon oscillation accuracies for 
spectroscopic redshift surveys as a function of survey area. We 
illustrate cases corresponding to the redshift windows that are 
naturally probed by optical spectrographs, 0.5 < z < 1.3 and 
2.5 < z < 3.5. We have assumed a sufficient number density of 
galaxies that shot noise is unimportant. A reference line Accu- 
racy (X Arca^-'^/'^ is plotted. The dependence of accuracy on area 
becomes steeper for small areas because the baryon oscillations 
are no longer being adequately resolved by the data. 



formation. Eisenstein et al. determined a 4% measurement 
of the acoustic scale when the clustering pattern was aver- 
aged over angles, using the full information contained in the 
shape. Our combined tangential and radial measurements 
suggest an overall accuracy of about 5% from just the oscil- 
latory component, which appears broadly consistent. 

5 CHANGING THE COSMOLOGICAL 
PARAMETERS 

These fitting formula coefficients have been derived from a 
grid of simulated surveys assuming a fiducial ACDM cosmol- 



ogy. However, the scaling arguments presented in Section 3 
apply more generally. As a result, it is a good approxima- 
tion to use the fitting formula of equation 6 for a range of 
cosmologicaJ parameters, if we compute the volume V, lin- 
ear growth factor D{z) and (where applicable) the radial 
position error ar using the new set of parameters. The co- 
efficients xo, Vq, (Jr,o and D{zq) should remain unaltered at 
their ACDM calibrations. However, two further changes are 
required: 

• The amplitude and shape of the input power spectrum 
P depend on the cosmological parameters. Our technique 
is largely insensitive to these dependences because we di- 
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Figure 8. The tangential baryon oscillation accuracy for spectro- 
scopic rcdshift surveys as a function of survey area for a series of 
redshift slices of width 0.2. We have assumed a sufficient number 
density of galaxies that shot noise is unimportant. A reference 
line Accuracy oc Area~^/^ is plotted. The dependence of accu- 
racy on area becomes steeper for small areas because tlic baryon 
oscillations are no longer being adequately resolved by the data. 
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Figure 10. The tangential baryon oscillation accuracy for pho- 
tometric redshift surveys as a function of survey area for a series 
of redshift slices of width 0.5. We have assumed a photometric 
redshift error ctq =0.03 and a sufficient number density of galax- 
ies that shot noise is unimportant. A reference line Accuracy oc 
Area"^/'^ is plotted. The dependence of accuracy on area be- 
comes steeper for small areas because the baryon oscillations are 
no longer being adequately resolved by the data. 
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Figure 9. The same as Figure 
cillation accuracy. 



vide out the overall power spectrum shape before fitting the 
baryon oscillations. However, the balance between cosmic 
variance and shot noise will be affected (i.e. the value of nP 
in equation 3). For a new set of parameters, the coefficient 
riefi should be scaled inversely with the characteristic power 
spectrum amplitude for the scales of interest, relative to its 
value in the fiducial case. 

• We should re-estimate the cut-off redshift Zm at which 
all of the high-amplitude acoustic peats become visible: the 
location of the non-linear transition scale at a given redshift 
depends on the growth of density perturbations, which is 
determined by the cosmological parameters. 



6 SPARSE-SAMPLING STRATEGIES 

Thus far, our formulae refer to surveys covering a fully con- 
tiguous sky area. However, the optimal strategy for measur- 
ing acoustic oscillations given a fixed observing time may 
not be to survey a contiguous area, but rather to sparsely- 



sample a larger area: gathering a larger density-of-states in 
Fourier space at the expense of an increased convolution of 
the input power spectrum (i.e. more smoothing of the acous- 
tic oscillations) and increased correlations between adjacent 
Fourier bins (i.e. less statistical significance for an observed 
peak or trough in power). In practice, sparse-sampling could 
be achieved by a non-contiguous pattern of telescope point- 
ing centres or, for a wide-field multi-object spectrograph, by 
distributing the fibres non-uniformly across the field-of-view. 

In the first approximation, the effectiveness of a sparse- 
sampling strategy depends on the angular size 9 of the ob- 
served survey patches (e.g. the field-of-view of the optical 
spectrograph) compared to the angular scale of the bary- 
onic features in the power spectrum ~ s/r(z) — 2.6 deg 
at a = 1. If ~ 1 deg then Wk will contain structure on 
scales similar to the acoustic preferred scale, and an unac- 
ceptable degree of convolution will result. If 9 -C 1 deg, then 
a sparse-sampling strategy will usually be preferred. 

We investigated this trade-off by simulating a series 
of sparse-sampling strategies with 9 = 2, 5, 10, 30 arcmin, 
considering spectroscopic redshift surveys only, and measur- 
ing the power spectr um in angle- a veraged bins of constant 
wavenumber k = \/k% -\- ky -\- k1. For each value of 9 we 
considered a series of survey 'filling factors' 1// such that 



/ = 



Sparsely sampled area 



(10) 



Observed area 

For the purposes of this simple investigation we assumed 
that the survey window function was a regular grid of square 
patches of size 9x9 (we note that other sampling strategies 
may be preferred, such as a random distribution of point- 
ings or a logarithmic spiral). For each (6',/) we determined 
an 'effective area gain' for the sparsely-sampled survey, by 
which we should multiply our observed (sparsely-sampled) 
area to produce the approximate input to the baryon oscil- 
lation fitting formulae. 

The Fourier transform grid required to analyze a vol- 
ume large enough to ensure a high-accuracy measurement of 
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Figure 11. EflFective area gain for a series of sparsely-sampled 
survey strategies at z = 1, varying the filling factor / of equation 
10 and the survey patch size d. The area gains are evaluated using 
equation 12 and are compared with the simple increase in survey 
performance neglecting the effects of convolution and mode cor- 
relations (the dashed line). The cases analyzed are indicated by 
the solid circles. 

the baryon oscillations, whilst maintaining a resolution sev- 
eral times better than the sparse sampling scale of a few ax- 
cmin, is prohibitively large. Therefore we adopted a different 
approach, estimating the effective area gain by quantifying 

three competing effects: 

• The average decrease in the amplitude A of the acoustic 
oscillations due to convolution with the window function. 

• The average decrease in the power spectrum error a in 
each Fourier bin (i.e. the diagonal elements of the covari- 
ance matrix) due to the increased number of Fourier modes 
analyzed. 

• The increased correlation of each Fourier bin with its 
neighbours, defined by quantifying an 'efiiective number of 
independent modes' vrn for each bin i using the covariance 
matrix dj of equation 1. For a uniform survey window func- 
tion in a cuboid, rrn = l/{dPif = {P + /Cu, where P is 
the power spectrum amplitude in bin i. For a general window 
function we defined: 



(11) 



such that the off-diagonal covariance matrix elements de- 
crease the independence of the bins. We take the sum over 
j up to the non-linear transition scale, and then define the 
average across the bins, m = rfii. 

We initially measured these quantities for a fiducial con- 
tiguous (/ = 1) survey of 100 deg^ spanning redshift range 
0.75 < z < 1.25. We then repeated our analysis for each pair 
of values of {9, f) defining in each case 



Effective area gain 



A/ a 
Ao/ao 



m 
mo 



(12) 



where the subscript '0' indicates values for the fiducial sur- 
vey. The relative powers of the quantities are chosen in accor- 
dance with their scaling with the number of Fourier modes 
m: Area oc m and a oc l/^/m. For the cases with small 
values of 6, the convolution involves small-scale power from 
the non-linear clustering regime, thus we modified our input 



linear power spectrum using the non-linear prescription of 
Peacock & Dodds (1994). 

The results are displayed in Figure 11. As expected, 
large survey patches ^ > 30 arcmin do not favour sparse- 
sampling strategies because of the consequent serious 
smoothing of the acoustic oscillations. If ^ < 10 arcmin then 
sparse-sampling strategies are preferred, although we note 
that the resulting performance plotted in Figure 11, which 
includes the window function effects, is not as good as 
that which would be inferred by using the entire 'sparsely- 
sampled area' as the input area in the fitting formula, thus 
neglecting the window function effects (as indicated by the 
'simple area gain' line plotted in Figure 11). We emphasize 
that our calculations here are only a first approximation and 
this is a subject requiring further study. 



7 SUMMARY 

We have developed a fitting formula for the accuracy with 
which the characteristic baryon oscillation scale may be ex- 
tracted from future spectroscopic and photometric redshift 
surveys in the tangential and radial directions, using heuris- 
tic scaling arguments calibrated using an accelerated ver- 
sion of the 'model-independent' method of Blalce & Glaze- 
brook (2003). The formula is given in equations 6 to 9 with 
the values of the parameters listed in Table 1, and repro- 
duces the simulation results with a fractional scatter of 7% 
(10%) in the tangential (radial) direction, over a wide grid of 
survey configurations. Simple modifications allow the fitting 
formula to be applied for a range of cosmological parame- 
ters. We have also investigated how a simple sparse-sampling 
strategy may be used to enhance the effective survey area 
if the sampling scale 9 is much smaller than the charac- 
teristic angular acoustic scale (6' ^ 1 deg). This may be 
implemented for a wide-field multi-object spectrograph by 
clustering the fibres in the field-of-view. 
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